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Abstract 

We compute the light hadron mass spectrum in quenched lattice QCD at (3 = 6.0 
using the Sheikholeslami-Wohlert fermionic action. The calculation is done for several 
choices of the coefficient csw, including csw = and the recently proposed optimal 
value csw = 1-769. We find that the individual masses change by up to 30% under 
0(a) improvement. The spectrum calculation suggests csw ~ 1-4 for the optimal value 
of the coefficient. 
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1 Introduction 



Wilson fermions have some advantages over staggered fermions. For example, one is not 
restricted to multiples of four flavors, the flavor symmetry is exact so that one does not have 
to test for its restoration [f], and the flavor and Dirac structures are not intertwined on the 
lattice making it easier to construct composite fermionic operators. Their disadvantage is 
that they are accompanied by cut-off effects of 0(a) , whereas the gauge field action and 
staggered fermions introduce corrections of 0(a 2 ) only. 

To be able to keep the size of the lattice at a reasonable level and yet explore large 
physical volumes, one is forced to do simulations at moderate values of the coupling, for 
example at (3 = 6/g 2 = 6.0 in the quenched theory. At this value the cut-off lies at O.f fm 
(taking the string tension as the physical scale), which is to be compared with the typical 
size of a hadron of f fm. The result is that one has to reckon with 0(a) effects causing 
systematic errors on the level of 20% or more. 

By adding the local counterterm 

l-Kgcsw(g)a^2^(x)cr llv F llv (x)^(x) (I) 

^ x 

(we assume r = I throughout the paper) to the Wilson fermionic action, cut-off effects in on- 
shell quantities [2] can be reduced from 0(a) to 0(a 2 ) [3]. Such quantities are, for example, 
hadron masses and energies as well as hadronic matrix elements of local operators. In the 
case of hadronic matrix elements the operators have to be improved as well. 

To achieve this, one has to compute the coefficient csw(g) to all orders in perturba- 
tion theory, that means non-perturbatively in practice. We call the fermionic action non- 
perturbatively 0(a) improved if csw is determined non-perturbatively, and we call the re- 
sulting value of csw an optimal value. Perturbatively csw is known to 0(g 2 ) [4]: 

csw(g) = 1 + 0.2659/. (2) 

The lowest order perturbative result, csw = 1, reduces cut-off effects to 0(a/lna) only, 
the second order result to (9(a/ln 2 a), and so on. A crude approximation of csw might be 
obtained by reorganizing the perturbative expansion such that the leading tadpole contri- 
butions are treated non-perturbatively [5], while the remainder is computed in low-order 
renormalized perturbation theory. This leads to 

csw(g) = u~ 3 (l + 0.0159/ 2 ), (3) 

where u = (|TrLVj} ? and g* is the coupling constant renormalized at some physical scale. 
Taking g* 2 = g 2 u^ 4 , we obtain at (3 = 6.0 the value csw = 1-518 for u = 0.8778. 

Recently, csw has been determined non-perturbatively in the quenched theory [6]. As 
a criterion for 0(a) improvement these authors use the PCAC relation, which connects the 
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divergence of the axial vector current with the quark mass and the pseudoscalar density, and 
demand that it be valid exactly so that every other on-shell quantity has possible corrections 
of 0(a 2 ) only. The result is [7] 

, , 1 - 0.656 g 2 -0.152 g A - 0.054/ , 

csw(g) = — n — , a < 1. 4 

yyj 1 - 0.922 # 2 ' y ~ y ' 

At small g this formula approaches the perturbative result (2). At (3 = 6.0 it gives csw = 
1.769, a value which is significantly higher than the tadpole improved perturbative result. 
There are also other possibilities to determine csw which we will address in a separate 
publication. 

The improvement program is a systematic expansion in powers of a, and one has to bear 
in mind that it is no longer useful if 0(a 2 ) corrections are of the same magnitude as 0(a) 
effects, which will occur as g increases. 

In this paper we shall look at the effect of improvement on the hadron mass spectrum 1 . 
Our first objective is to see how sensitive masses and mass ratios are to the values of csw- 
Secondly, we would like to investigate whether a determination of csw from the mass spec- 
trum is possible. 



2 The Calculation 

This work is the first step in the aim to apply 0(a) improvement to our structure function 
calculations [9]. 

We work at (3 = 6.0 and use lattice sizes 16 3 32 and 24 3 32. In most of our runs we 
generated 100 - 200 configurations. For the pure Wilson case we have accumulated (9(5000) 
configurations at k = 0.155, 0.153 and 0.1515 on the 16 3 32 lattice in our attempt to determine 
the gluon distribution functions of the nucleon [10]. For the gauge field update we use a 
combination of 16 overrelaxation sweeps followed by a three-hit Metropolis update. This 
procedure is repeated 50 times to generate a new configuration. 

The implementation of the improvement term (1) in our calculation is described in 
ref. [11]. For the matrix inversion we mainly used the minimal residue algorithm, except for 
the lightest quark mass on the larger lattice where we used the BiCGstab algorithm [12]. 
At csw = 3, our largest csw value, we experienced convergence problems with the minimal 
residue algorithm. So we experimented with the BiCGstab inversion algorithm, but this did 
not seem to improve the situation significantly. 

We used Jacobi smearing for source and sink with 50 iterations and a smearing parameter 
k s = 0.21 [13], giving a radius of the hadrons of about four lattice spacings. We employed 

1 For a similar investigation at (3 = 5.7, which reached us after this paper was completed, see ref. [8]. 
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c S w = 


V 


K 


m v 


nip 




ao 


m ai 


m bi 


16 3 32 


0.1487 
0.1515 
0.1530 
0.1550 


0.637(2) 
0.5033(4) 
0.4221(4) 
0.2966(5) 


0.682(2) 
0.5682(7) 
0.5058(8) 
0.4227(15) 


1.061(9) 
0.902(2) 
0.798(2) 
0.652(3) 


0.891(24) 
0.817(7) 
0.763(11) 
0.735(15) 


0.940(21) 
0.851(7) 
0.797(6) 
0.717(12) 


0.943(32) 
0.849(13) 
0.809(7) 
0.736(9) 


24 3 32 


0.1550 
0.1558 
0.1563 


0.292(2) 
0.229(2) 
0.179(3) 


0.418(5) 
0.384(7) 
0.358(11) 


0.638(8) 

0.555(12) 

0.488(22) 


0.610(48) 
0.616(90) 
0.88(15) 


0.657(33) 
0.613(41) 
0.584(52) 


0.659(35) 
0.638(38) 
0.615(44) 


c. 1. 


0.15717(3) 





0.328(5) 


0.421(15) 


0.656(19) 


0.622(15) 


0.646(14) 



Table 1: The hadron masses in lattice units for pure Wilson fermions with csw = 0. In the 
bottom row we give n c and the mass values extrapolated to the chiral limit. The numbers 
in roman (italic) are from three-parameter (two-parameter) fits. 









csw 


= 1.769 








V 




m v 


m p 


m N 


m ao 


m ai 




16 3 32 


0.1300 
0.1310 
0.1320 
0.1324 
0.1333 
0.1342 


0.701(4) 
0.624(4) 
0.545(5) 
0.501(2) 
0.410(3) 
0.299(3) 


0.782(5) 
0.715(5) 
0.644(8) 
0.613(4) 
0.547(5) 
0.488(10) 


1.17(2) 

1.062(17) 

0.974(16) 

0.906(13) 

0.791(16) 

0.674(27) 


0.785(17) 

0.74(3) 

0.85(8) 


0.817(12) 
0.750(15) 
0.683(20) 


0.831(23) 
0.778(17) 
0.744(26) 


24 3 32 


0.1342 
0.1346 
0.1348 


0.302(2) 
0.238(2) 
0.188(7) 


0.492(5) 
0.470(9) 
0.457(20) 


0.695(9) 

0.642(16) 

0.605(22) 


0.82(3) 
1.00(8) 
1.52(20) 


0.715(19) 
0.684(26) 
0.664(34) 


0.758(16) 
0.745(20) 
0.736(29) 


c. /. 


0.13529(3) 





0.429(11) 


0.548(25) 


0.845(38) 


0.629(19) 


0.716(18) 



Table 2: The hadron masses in lattice units for non-perturbatively 0(a) improved Wilson 
fermions with csw = 1-769, according to eq. (4). In the bottom row we give n c and the 
mass values extrapolated to the chiral limit. The numbers in roman (italic) are from three- 
parameter (two-parameter) fits. 
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c sw = 1-92 


V 


K 


m v 


m p 


m N 


m ao 


m ai 


m bl 


16 3 32 


0.1290 
0.1300 
0.1310 
0.1320 


0.644(1) 
0.553(1) 
0.451(2) 
0.328(2) 


0.734(2) 
0.661(3) 
0.586(4) 
0.516(10) 


1.121(9) 
1.001(10) 
0.875(11) 
0.735(20) 


0.888(9) 
0.819(11) 
0.755(14) 
0.793(65) 


0.946(10) 
0.879(11) 
0.812(13) 
0.763(29) 


0.959(11) 
0.892(12) 
0.827(14) 
0.771(34) 


c. 1. 


0.13322(6) 





0.430(32) 


0.539(88) 


0.611(33) 


0.668(27) 


0.681(30) 



Table 3: The hadron masses in lattice units for improved Wilson fermions with csw = 1-92. 
In the bottom row we give n c and the mass values extrapolated to the chiral limit. The 
numbers in roman [italic) are from three-parameter (two-parameter) fits. 



csw = 2.25 


V 


K 


m v 


m p 


mjv 


m ao 


m ai 


m bl 




0.1250 


0.639(1) 


0.739(3) 


1.097(10) 


0.885(16) 


0.940(22) 


0.946(24) 




0.1260 


0.535(2) 


0.658(4) 


0.966(12) 


0.800(18) 


0.848(22) 


0.861(25) 


16 3 32 


0.1265 


0.476(2) 


0.617(4) 


0.893(11) 


0.774(26) 


0.802(26) 


0.808(33) 




0.1270 


0.414(2) 


0.576(5) 


0.807(15) 


0.738(38) 


0.754(26) 


0.750(43) 




0.1277 


0.305(3) 


0.514(12) 


0.684(19) 




0.733(61) 


0.697(33) 


c. 1. 


0.12861(5) 





0.435(32) 


0.461(80) 


0.597(56) 


0.602(49) 


0.597(44) 



Table 4: The hadron masses in lattice units for improved Wilson fermions with csw = 2.25. 
In the bottom row we give n c and the mass values extrapolated to the chiral limit. The 
numbers in roman [italic) are from three-parameter (two-parameter) fits. 



csw — 3 


V 


K 


m v 


m p 


m N 




0.1150 


0.705(4) 


0.840(4) 


1.239(9) 




0.1155 


0.643(4) 


0.783(4) 


1.172(9) 


16 3 32 


0.1160 


0.555(6) 


0.727(3) 


1.054(9) 




0.1165 


0.470(12) 


0.684(6) 


0.936(13) 




0.1170 


0.408(7) 


0.603(7) 


0.858(12) 




0.1173 


0.342(12) 


0.569(13) 


0.758(16) 


c. /. 


0.11826(13) 





0J f 64(13) 


0.558(25) 



Table 5: The hadron masses in lattice units for improved Wilson fermions with csw = 3. In 
the bottom row we give n c and the mass values extrapolated to the chiral limit. The latter 
values are from two-parameter fits. 
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K c 






eci. (5) 

1 v / 


eq. (Q) 


eq. (7) 





0.15693(1) 


0.15734(5) 


0.15717(3) 


1.769 


0.13522(1) 


0.13534(4) 


0.13529(3) 


1.92 


0.13306(2) 


0.13333(11) 


0.13322(6) 


2.25 


0.12850(2) 


0.12866(8) 


0.12861(5) 


3 


0.11799(4) 


0.11862(39) 


0.11826(13) 



Table 6: The critical values of k, k C} for the linear (eq. (5)), chiral (eq. (6)) and phenomeno- 
logical fit (eq. (7)) for our various csw parameters. 



both relativistic and non-relativistic wave functions [9, 13], except for the high statistics runs 
where we only looked at the non-relativistic wave function in order to save computer time. 

We calculated the masses of 7r, /?, nucleon, a , a\ and b\. In the following we shall mainly 
focus on the 7r, p and nucleon masses. The csw values we considered were (pure Wilson), 
1.769 (the value of eq. (4)), 1.92, 2.25 and 3. In addition, we shall make use of the UKQCD 
results [14] for the ir and p mass at csw = 1.4785, which is their estimate of the tadpole 
improved perturbative value. 



3 Results 

We present our results in tables 1-5. The errors on the data are bootstrap errors. For the 
meson masses we found very little difference between using relativistic and non-relativistic 
wave functions, and we settled for relativistic wave functions (except for the high statistics 
runs). For the nucleon we have chosen non-relativistic wave functions which performed 
slightly better. 

To obtain the critical value of k, k C} and the hadron masses in the chiral limit, we 
extrapolate our data to zero ir mass. We expect 

ml = b(^-^). (5) 

Using this relation gives a rather poor fit of the data, and we can see that there is slight 
curvature in a plot of m\ against 1/k. Quenched chiral perturbation theory predicts [15] 

m ; = *(!_!)*, ( 6 ) 
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0.16 




0.1 1 

0.0 1.0 2.0 3.0 

c sw 



Figure 1: The critical value of k against csw- The points at csw = and 1.769 (•) include 
data at small ir masses on the 24 3 32 lattice, the points at 1.92, 2.25 and 3 (□) are from the 
16 3 32 lattice. The data point at csw = 1.4785 is from UKQCD [14]. The solid curve is the 
prediction of perturbation theory, the dashed curve that of tadpole improved perturbation 
theory. The dotted line connecting the data points is meant to guide the eye. 



7 




8 



1 .6 



m /m 
N p 



1 .5 



1 .4 



1 .3 



1 .2 




0.0 



0.2 



0.4 



0.6 



0.8 

(m /m )2 
71 p 



1.0 



Figure 3: APE plot for csvf = (A) and csvf = 1.769 (•). The errors are bootstrap errors. 
The solid boxes represent the experimental value and the expected result in the heavy quark 
limit, respectively. 
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where 8 is small and positive. We made fits using this formula but found that 8 was always 
negative. We concluded that our k values are too far from n c for the formula to be applica- 
ble. This is in agreement with observations made by other authors [16]. As an alternative 
parameterization of the curvature we used the phenomenological fit 

- = — + b 2 ml + b 3 m 3 v . (7) 

In table 6 we give the values of n c for the different fits. The linear fits give \ 2 /d.o.f. values 
of up to 20. The other two fits give both acceptable values of % 2 , but the phenomenological 
fits give slightly better results. In the following we shall take n c from the phenomenological 
fits. 

In fig. 1 we plot k c against csw- We observe a rather linear behavior of n c with csw- It 
is interesting to compare the numerical values with the perturbative result 

k c = -[1 + g 2 (0. 108571 - 0.028989 c sw - 0.012064 c| w )] (8) 

8 

and tadpole improved perturbative result 

u k c = hi + g* 2 (0. 025238 - 0.028989 c sw u 3 - 0.012064 (c S wu 3 ) 2 )]. (9) 
8 

Neither formula agrees well with the data, but the tadpole improved line gives rough agree- 
ment in the neighborhood of csw = 1-769. 

We fit the rest of the hadron masses by the formula 

m 2 H = m° H 2 + b 2 ml + 6 3 m^. (10) 

This gives a better fit to the data than the ansatz [17] 

m H = m° H + b' 2 ml + b' 3 ml. (11) 

Note that the two formulae differ only by terms of 0(m^). For csw = 3 as well as for the 
a , ci\ and b\ masses only a two-parameter fit with b 3 set to zero was reasonable. The mass 
values in the chiral limit are also given in tables 1-5. We consider our extrapolations 
for csw = and 1.769 to be most reliable, because on the larger volume we could explore 
relatively small ir masses. 

We show in fig. 2 the 7r, p and nucleon masses for csw = and 1.769 together with the 
fits. It is important to have many k values in the fit to do a reliable extrapolation. At the 
third smallest quark mass we have results on two different volumes. The values agree within 
the errors. This indicates that all our results on the 16 3 32 lattice do not suffer from finite 
size effects. 

In fig. 3 we show the APE plot for csw = and 1.769 together with our chiral extrapo- 
lations. We find rather different results for the two csw values. However, in the chiral limit 
the nucleon to p mass ratio is about the same. The csw = 1-769 result appears to go more 
smoothly to the heavy quark limit. 
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1.6 



m / m 
N p 



1 A 



1 .2 



1 .0 



0.8 



0.0 



1 .0 



2.0 3.0 

c sw 



Figure 4: The nucleoli to p mass ratio in the chiral limit against csw- The symbols are as 
in fig. 1. The dashed line is the experimental value. 
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0.50 



m 
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0.45 



0.40 



0.35 



0.30 



0.0 



1 .0 



2.0 



3.0 



sw 



Figure 5: The p mass in the chiral limit against csw- The symbols are as in fig. 1. The 
dashed line is the experimental value, where the lattice spacing is taken from the string 
tension: K = 0.0515(28) [18], with K = (427 MeV) 2 as given by the Cornell potential [19]. 
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As we noted in the APE plot, the nucleon to p mass ratio in the chiral limit was the 
same at csw = and 1.769. In fig. 4 we show the ratio for all our csw values. The results 
agree with the experimental value within the errors and seem to be independent of csw- 
Away from the chiral limit the mass ratio depends more strongly on csw- This means that 
if one is only interested in mass ratios in the chiral limit, then Wilson fermions seem to give 
reasonable extrapolations. 

Although the nucleon to p mass ratio is independent of csw, when we look at the p 
(and nucleon) mass itself in the chiral limit we see a strong csw dependence, as is shown in 
fig. 5. Varying csw from to 1.769 we observe a 30% increase of the mass. Even between 
the tadpole improved perturbative result csw ~ 1-5 and csw = 1-769 the mass changes by 
about 10%. Expressing the p mass in terms of the string tension K, which has discretization 
errors of 0(a 2 ) only, gives the dashed line. If one assumes that the effect of quenching can 
be neglected, one can use these results to estimate the optimal value of csw up to possible 
0(a 2 ) corrections. 

4 Discussion 

We have done a systematic investigation of the csw dependence of the light hadron mass 
spectrum. We see that at (3 = 6.0 it is important to choose a good value of csw- 

Disregarding errors due to quenching and 0(a 2 ) corrections, our best guess for the optimal 
value of csw is csw ~ 1-4. This result is somewhat lower than the value given by eq. (4). If 
the difference is attributed to 0(a 2 ) effects, we would estimate that 0(a 2 ) corrections can be 
half as big as 0(a) corrections at (3 = 6.0. This would mean that the on-shell improvement 
program does not make much sense at much lower values of (3, for example at (3 = 5.7. 

We consider our work only to be a first step. Computing masses in the chiral limit 
is difficult. A lesson we have learned is that one needs many k values for a reliable chiral 
extrapolation. A more precise determination of the optimal value of csw would require much 
higher statistics. 
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